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This  Technical  Report  was  prepared  by  the  Mechanical  Behavior  and  Aging 
Section  (MKPB)  of  the  Air  Force  Rocket  Propulsion  Laboratory  under  Job 
Order  No.  2307M1KT.  The  principal  author  was  Michael  H.  Quinlan,  now  at 
the  University  College  Cork,  Cork  City,  Ireland.  Some  sections  of  the 
report  (as  detailed  in  the  Background  Section)  were  written  by  Durwood  I. 
Thrasher. 


This  report  has  been  reviewed  by  the  Information  Office/XOJ  and  is  releasable 
to  the  National  Technical  Information  Service  (NTIS).  At  NTIS  it  will  be 
available  to  the  general  public,  Including  foreign  nations.  This  report 
is  unclassified  and  suitable  for  general  public  release. 
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'whenever  c = (<j>  ( 1 -ct)  + ay)  / au  = 0,  and 
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1.0  BACKGROUND 


This  report  is  the  interim  product  of  more  than  eight  years  of  research 
by  Dr.  Michael  Quinlan.  He  was  sponsored  and  provided  with  facilities  by 
the  AFRPL  during  the  period  October  1975  through  October  1977  as  a National 
Research  Council  Postdoctoral  Research  Associate.  From  October  1977  through 
April  1979,  Dr.  Quinlan  continued  work  on  the  theoretical  development  and 
analysis  of  propellant  test  data  (generated  during  his  work  at  the  AFRPL) 
at  University  College  Cork,  Ireland. 

Equipment  and  propellant  samples,  and  some  technical  support  for 
Dr.  Quinlan's  experimental  efforts,  were  provided  under  AFRPL  In-House 
Project  57301 3NE,  Solid  Propellant  Mechanical  Behavior.  Subsequent  support 
was  provided  under  AFRPL  Project  230/M1KT,  Propellant  Constitutive  Theory, 
until  Dr.  Quinlan  returned  to  Ireland  in  1977.  His  work  since  that  time  was 
conducted  without  support  from  the  Air  Force  by  mutual  agreement. 

Dr.  Quinlan  authored  Section  3,  Theoretical  Development.  Section  3 
comprises  the  major  content  of  this  report.  The  remainder  of  the  report 
was  authored  by  the  current  Air  Force  Project  Manager,  Durwood  Thrasher. 
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2.0  tXPERIMENTAT ION 


Varied  experiments  were  performed  by  Dr.  Quinlan  during  his  research 
at  the  AFRPL  on  two  propellants:  ANB-3066  (a  carboxy- terminated  polybutadiene 
(CTPB)  propellant  containing  ammonium  perchlorate  and  aluminum)  and  TP-H1011 
(a  polybutadiene— acryl  ic  acid-acrylonitri le  (PBAN)  propellant  containing 
ammonium  perchlorate  and  aluminum). 

Much  of  the  testing  was  done  to  gain  qualitative  results  on  the  general 
mechanical  behavior  of  the  two  solid  propellants  as  well  as  to  develop  a 
characterization  technique  to  allow  determination  of  a specific  nonlinear 
constitutive  law. 

The  test  specimen  used  is  portrayed  in  Figure  1(a).  The  specimen  is 
a cylindrical  uniaxial  tensile  bar,  carefully  milled  to  achieve  a uniform 
diameter.  The  ends  of  the  specimen  were  cut  using  a razor  blade  and  a 
special  jig  to  obtain  square  ends  and  a reproducible  length.  The  ends  of 
the  specimen  were  bonded  to  aluminum  end  tabs  using  epoxy  adhesive.  A jig 
was  used  to  assure  accurate  alignment.  The  aluminum  end  tabs  were  designed 
to  fit  into  special  grips  and  incorporated  set  screws  to  clamp  the  end  tabs 
secu-ely.  This  feature  prevented  slipping  of  the  specimen  in  the  grips  under 
load  reversal.  Considerable  effort  was  expended  to  eliminate  play  and  lost 
motion  in  the  load  linkage  of  the  test  machine.  Some  troublesome  "strange" 
,'pellant  behavior  was,  in  fact,  traced  to  elastic  "play"  in  a part  of  the 
load  linkage. 

The  testing  was  performed  using  a standard  laboratory  tensile  tester, 
manufactured  by  Instron  Corporation.  The  test  assembly  included  a gas 
dilatometer  manufactured  by  Dr.  Richard  J.  Farris,  and  the  dilatometer  load 
cell  was  used  to  measure  force  on  the  specimen.  The  specimen  extension  was 
determined  from  crosshead  travel  of  the  tensile  tester.  A temperature 
conditioning  unit  was  used  to  maintain  a stable  "room  temperature"  77°F 
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environment  on  the  specimen.  (While  problems  with  the  dilatometer  precluded 
dilatation  measurements,  its  load  linkage  was  found  to  be  very  satisfactory 
for  obtaining  data  under  strain  histories  that  produced  load  reversals.) 

The  characterization  test  eventually  selected  is  schematically  shown  in 
Figure  1(b).  The  test  consists  of  a series  of  sawtooth  events  of  increasing 
strain  amplitude.  A long  (several  minutes)  "rest"  period  at  zero  strain 
is  interposed  between  the  sawtooth  events  to  allow  the  visoelastic  part 
of  the  propellant  response  to  decay  to  a negligible  level.  The  loading  and 
unloading  parts  of  each  sawtooth  have  the  same  constant  strain  rate. 

Some  useful  variations  of  the  above  strain  history  include: 

a.  Insertion  of  a sawtooth  event  with  a lower  maximum  strain 


than  the  previous  maximum. 

b.  Insertion  of  a relaxation  process  between  sawtooth  events. 

c.  Changing  the  strain  rate  for  specific  sawtooth  events. 

Data  from  the  characterization  test  (as  described  above)  performed  on 
a single  TP-H1011  propellant  sample  were  used  to  produce  the  propellant 
constitutive  law  constants  described  at  the  end  of  the  following  section 
of  the  report.  The  overall  response  of  the  TP- HI 01 1 specimen  was  similar 
to  that  of  all  TP-H1011  specimens  tested.  The  ANB-3066  specimens  also 
showed  the  same  general  response  characteristics  (e.g.,  rejoining  the 
virgin  stress-strain  curve  when  strained  to  a level  higher  than  the 
previous  maximum  strain).  These  characteristics  have,  of  course,  been 
widely  observed  and  reported  for  solid  propellants  in  general. 
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3.0  DEVELOPMENT  AND  APPLICATION  OF  THEORY 


3 . 1 Introduction 

We  consider  the  behavior  of  cylindrical  solid  propellant  bodies, 
of  small  diameter- to-height  ratio  and  weight,  subject  to  small  deforma- 
tions which  preserve,  and  are  parallel  to,  the  axis.  A temporal 
sequence  of  deformations  emanating  from  the  same  equilibrium  state  at 
time  zero  and  terminating  at  some  individual  current  time  is  called  a 
process  and  the  average  traction  subsiding  on  the  disk  boundaries  upon 
termination  is  called  the  force-response.  The  value  of  a deformation 
in  a process  at  a given  time  is  the  axis  length  at  that  time  divided 
by  its  initial  length.  Assuming  that  the  force-response  to  any  given 
process  is  unique,  our  task  is,  then,  to  adequately  describe,  or  model, 
this  correspondence. 

Structurally,  propellants  consist  of  minute  (i.e.,  in  the  order 
of  microns)  rigid  particles  densely  embedded  in  a rubbery-like  polymer. 
Such  materials  respond  to  processes  wi  th  an  alteration  in  the  amount  of 
species  to  species  bonding  - hence  the  idea  of  a process  inducing  or 
creating  a particular  current  state  of  bonding.  An  adequate  continuum 
description  of  bonding  changes  can  be  effected  by  the  introduction  of  an 
equivalance  relation  on  the  totality  of  processes.  The  intuitive 
interpretation  attached  to  two  processes  being  equivalent  is  that  they 
both  induce  the  same  state  of  species  to  species  bonding  at  their 
respective  current  times.  The  validity  of  taking  such  a concept  as  being 
basic  can  only  be  measured  by  our  ability  to  frame  intuitively  plausible 
and  useful  axioms;  useful  in  the  sense  of  ultimately  unravelling  the 
force- response  function.  Of  course  if  one  is  persuaded  to  accept 


(6,1)  and  (6,2)  as  being  basic  on  the  basis  of  matching  experimental  data, 
then  all  considerations  related  to  bonding  can  be  ignored.  This  would  be 
as  unjustified  as  to  accept  the  Navier  equation  of  elasticity,  with  all 
boundary  conditions  in  terms  of  displacements,  and  ignore  the  concept  of 
stress. 

Allied  to  the  notion  of  a state  of  bonding  is  the  idea  of  continuing 
or  extending  a process  for  any  preassigned  length  of  time  so  that  the 
extended  and  original  process  create  the  same  state  of  bonding.  As  the 
preassigned  lengths  of  time  are  allowed  to  increase  without  bound,  we 
suppose  that  the  corresponding  current  stretches  tend  to  a limit.  This 
limit,  called  a bending  pciumcdeA,  serves  to  label  the  state  of  bonding 
created  by  the  original  process. 

The  physical  requirement  that  a unique  state  of  bonding  subsist  in 
each  equilibrium  state  is  equivalent  to  the  requirement  that  all  processes 
with  the  same  bonding  parameter  create  the  s ante  state  of  bonding  and 
conversely.  Assuming  that  such  is  the  case,  the  set  of  all  processes 
with  the  same  bonding  parameter,  and  hence  each  state  of  bonding,  can  be 
identified  by  a unique  number,  namely,  the  bonding  parameter  of  any,  and 
hence  all,  processes  in  the  set.  Correspondences  of  this  type  are  not 
ji, common  in  mechanics,  indeed,  material  particles  are  defined  by  the 
requirement  that  they  can  be  put  into  one-one  correspondence  with  points 
in  a Euclidean  space.  Further  associations,  while  important,  are  entirely 
subjecti ve. 

These  ideas  suitably  axiomatized,  serve  as  a basic  framework  for 
acquiring  an  explicit  knowledge  of  the  evolution  of  bonding  states 
(i.e.,  the  bonding  parameter  at  each  time  versus  time)  in  a process. 


Motivated  to  some  extent  by  reaction  rate  theory  and  the  fact  that  various 
integral  representations  are  inadequate,  we  explore  the  possibility  of 
representing  the  various  evolutions  of  bonding  states  as  the  corresponding 
solutions  of  an  ordinary  differential  equation. 
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3.2  Basic  Structure 


Let  an  open  neiohbourhood  D of  unity  in  the  strictly  positive 

reals  be  qiven.  Then  a process  <p  of  duration  d > 0 is  a continuous 

ip 

function 

P : [O.d^J  - D ; p(0)  = 1 

whose  right  derivative  ^ is  piecewise  continuous.  An  extension  p(°) 
of  a process  p is  called  a o-  static  continuation  provided  that 
(<0 

ip  (t)  = p(d^)  for  all  t £ (d^.d^+o)  where  d^+o  is  the 

duration  of  p^0)  The  notation  p^  for  t € [0,d^)  is  used  to  denote 
the  process  gotten  by  restricting  p to  [ 0 , t ) . 

A process  p gotten  by  extending  a process  p such  that  p^  ~ ^ for^ 
all  t £ ( d^ , d^ ) is  called  a frozen  bonding  continuation  of  ip.  On 
the  question  of  whether  each  process  can  be  so-continued,  and  if  so, 
for  how  long?,  we  assert  that 

AM.,  Each  process  p has  a unique,  frozen  bonding  continuation  p5  of 
duration  d +6  for  each  <5  > 0 and  that  the  limit 

n(ip)  d=f  Lim^  cp5  ( d^+6 ) 

exists. 

We  note  that  AM?  is  a strengthened  version  of  axiom  3 in  [II. 

I he  statement  that,  a ciiq,.  state  of  bonding  subsists  in  a given 
equilibrium  configuration  means  that  the  (generally  different)  frozen 
bonding  continuations  of  ip  and  p seek  the  same  ultimate  stretch 
whenever  p ~ p and  conversely.  Mathematically  this  statement  is: 

AM , ( b ) i et  p and  p be  processes.  Then 

<p  — p i ft  'i(.p)  - ir(<p) . 

Clearly,  a given  state  of  bonding  is  induced  or  created  by  some 
process  p and  all  processes  tp  that  create  the  same  bonding  as  p 

l<p  ~ p means  that  ip  and  p create  the  same  state  of  bonding. 
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satisfy  n(ep)  n (m>)  n in  which  case  n ilur.ti  tn  i;i>s  t h«»  t t «' 
uniquely,  The  positive  number  n(cp)  is  called  the  \ j ts -* / < »• 

corresponding  to  the  process  cp.  It  Co  Hows  that  axiom  part  (b) 
in  [1]  is  an  immediate  consequence  of  AM and  AMj(b). 


Clearly,  the  bond i no  state  created  by  ' ouoht  to  become 
constant  as  o increases  without  bound  in  which  case  ^ is 

essentially  a frozen  bonding  continuation  of  g/'  for  all  o 
sufficiently  large  and  each  a > 0. 

AM j (n ) Let  cp  be  any  process,  lhen 

) ■*  g>(d  ) smoothly  as  o 


The  quail fication  smoothly  is  token  to  mean  that  the  derivative  of 
o m ) exists  and  limits  to  zero  as  o -♦  «■.  A knowledge  of 

events  prior  to  time  zero  is  irrelevant  iff  an  equilibrium  state 
subsists  at  time  zero,  let  s (cp)  be  the  force-response  (stress)  to 
the  process  cp  and  lot  the  statement  that  cp  is  a r:'.;;;f  f isoie/.sf  i\"i 
of  the  process  cp  be  taken  to  mean  that  <p(t)  - >p(tfa)  for  all  (Mvd^ 
and  (p  *>  1 on  (0,a)  and  cU"  d^f  a for  some'  eWi.  lhen  the  assumed  exist- 
ence of  an  equilibrium  state  at  time  zero  is  entirely  equivalent  to: 

£ lei  tp  ho  a right  translation  of  the  process  >t\  lhen  vp  vp  and 
s(>;’)  ( pi 

We  note  that  In  11],  t is  implicit  in  the  definition  of  a process  class. 


I ot  the  processes  cp  and  cp  have  identical  duration  times  lhen 
the  (pointwiso)  difference  cp  - cp  and  the'  distance 
d(cp,i|>)  supl  | (t)  «,'( ♦ ) ( tclO.d^))  an’  writ  dc’fine.l.  It,  sa\  , 

d*  - cl  we  form  the'  right  translation  r of  cp  with  K'  d , and  take' 

the'  distance  from  cp  to  cp  to  be  d(cp.cp) . tilearly  the'  totality  o*  pis' 
cesses  has  t tic'  structure  of  a psuodo  metric  space  with  respect  to  d. 
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It  is  simple  matter  to  recast  the  theory  of  materials  with 
variable  bonding  as  presented  in  Ml  with  axiom  1 ( resp . 2 ) replaced 
by  the  weaker  (resp.  stronger)  AM j (resp.AMp) , and  verify  that  the 
various  definitions,  assumptions,  and  propositions,  if  not  un- 
altered, have  obvious  content  preserving  modifications. 


3.3  Bond  breakage  and  I urination  as  a Rate  Process 

The  poasibiliiy  of  representing  bond  breakage  and  formation  as 
a rate  process  is  mathematical ly  equivalent  to  finding  a function  P 
such  that  the  bonding 

t C (0,d^)  - n(g>t)  (3,1) 

corresponding  to  any  process  ip  may  be  characterized  as  the  unique 
function  satisfying 

= P(n(f>t)*  <P(t),  ‘b(t))  v 0 < t < d (3,?) 

together  with  the  limit  condition 
n(tj^ ) ■♦  1 as  t t 0 

by  virtue  ot  i.  Wo  tacitly  assume  that  (3,1)  has  the  smoothness 
properties  ot  a process  on  (0,d  ).  Clearly,  then,  (3,1)  has  a unique 
continuous  extension  to  [0,d  ) and  the  notation  '»(m’0)  is  meaningful. 
Ihroughout  the  remainder  of  this  section  we  let  4)  be  a generic  process. 

Subject  to  an  a fortiori  verification,  we  proceed  upon  the 
supposition  'hat  such  .1  continuous  function  !’  exists.  ( I wly  AM. 
p.u  t (a)  and  the  continuity  of  P ensure  that  l’(n(gi)  v(<^),0)  =■  0. 
Conversely,  suppose  that  the  value  of  P at  (n(>r),  c(d  ),<'}  iS  zero, 

(p  ] . 

then  s'  ,uni  the  constant  funct  ion 


1 j 


(ft 

> 


r 
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i t [d  ,d  4o)  -»  n(u>) 

4'  V 

satisfy  (3,?)  on  |d  ,d  4o)  for  each  a > 0.  Since  this  solution  nu‘>t 
be  unique  and  AM^  part  (a)  must  be  satisfied,  wr  infer  that 
n(ip)  ip(d^).  Clearly,  then,  the  restriction 

P ( n (tp)  • »p(d  ),0)  = 0 iff  n(<p)  <p(d  ) (3.3) 

ip  »p 

must  always  be  satisfied. 

Since  each  process  must  have  a frozen  bonding  continuation,  it 
follows  from  (3,2)  that  the  algebraic  equation 

P(n(tp),  ip(  dp)  • • ) - 0 

must  have  at  least  one  root  o(tp).  We  require  that  r(tp),  called  the 
neutral  rate  corresponding  to  vp,  ba  unique. 

For  tp  to  have  a unique  frozen  bonding  continuation  up  to  time 
d^+fl,  it  is  both  necessary  and  sufficient  that  a unique  function 

t.  e (d^,d  4,s)  z(t)  exists  such  that 

P(n(,p),  Z (t).  >(t))  = 0;  z(d^)  = <p(dj;  z(d^)  - e(ip)  (3.4) 

for  all  t € Id  ,d  4,1),  While  necessary  renditions  that  ensure  exist  erne 
IP 

ire  not  < urrenf ly  available,  a plausible  sufficiency  condition,  that 
ensures  existence  for  all  sufficiently  small  S > 0,  is,  namely 

P(n(ip),  io(d^),-)  has  a non  zero  derivative  at  e(v’)  (3 ,S) 

this  condition,  demandinu  that  tlie  rate  of  change  of  bonding  be  locally 
directed  by  the  process  rate,  suffices  to  establish  the  existence  of 
on  oju-n  neighbourhood  0 of  m»(  »t  and  unique  continuously  ui ffeivnt tab  !r 
l unct ion 

i ('<  («p) , • ) ; o c 0 -»  R 
called  the  neutral  such  that 
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(a)  C(n(u>) , *P(cy ) - e(<p) 

(b) each  forward  solution  to 

i - C(n(g>).z)  ; z(d^)  * ip(d^)  (3,6) 

satisfies  (3,4)  uniquely. 

We  require  that  (3,5)  be  satisfied,  thereby  ensuring  that  >p  has  a 
unique  frozen  bonding  continuation  up  to  d^a  for  all  s sufficiently 
small.  It  follows  that  AM?  is  satisfied  if  (3,6)  has  a solution 
on  It,-)  (cf  sec  3.  4). 

A Taylor  series  expansion  of  P(n(g>t),  g>(t),-)  about  P(^) 
yields  the  first  order  approximation: 

ft(g>t)  ■ a(n((pt),  q>(t))  (iJ>(t)-C(n(ipt) , <p(t))]  ^,7) 

In  view  of  (3,3) 

^(n('Pj)»  s*>( t ) ) ■ 0 iff  tp(t)  » n(ip^)t  (3.8) 

We  assume  that  a is  the  constant  function  with  value  a. 

To  avoid  a return  to  the  classical  theory  of  viscoelasticity,  modulo 
a driving  term  ^ (c.f.[l],  sec.  9)  we  assume  that  a / 1.  Finally  in 
view  of 

*(>Pt)  - « J»(t)  when  n(»pt)  * »p(t) 
we  impose  the  growth  conditions  f)  < < 1 . 


-14- 


3.4  Frozen  Bonding  Continuations;  The  Neutral  Koto  I unit  ion 

Let  the  processes  xp  and  xp  have  the  same  current  bonding 
parameter  nQ  and  satisfy 

0 < nQ  = «P(d^)-"0  = (*M) 

For  significant  values  of  n0*  experience  of  bond  reformation  (resp. 
rupture)  indicates  that  each  frozen  bonding  continuation  of  ip 
(resp.  xp)  ought  to  be  weakly  increasing  (resp.  significantly 
decreasing)  in  each  sufficiently  small  right  neighbourhood  of  d . 

This  implies  that  p(xp)  = c(n0»~n0+n0)  (resp.  e(xp)  = C(n0»'\,+"0)) 
ought  to  be  small  and  positive  (resp.  significantly  negative). 

Clearly,  then,  the  function 

n -*  C(nQ,n+no)  d§f  C(nQ,n) 

is  defined  in  a neighbourhood  N of  zero  and  is  plausibly  of 
exponential  character  in  the  sense  that: 

NRF  There  exists  an  upper  bound  n(no)  such  that  the  derivative  of 
n c(  A0<  n)-  p(nQ)  is  proportional  to  itself,  the  constant  of 
proportional i ty  u(n0)  being  strictly  positive. 

Immediately,  then, 

C(«0.n)  -p(nQ)  - c(no)  v n i N. 


As  •>  consequence  of  (3,8),  it.  follows  that  c(n  f)  ,i(no). 

Moreover  if  u(n^) , p("0)  and  o(n^),  p(n^)  satisfy  NRF  , then  the 
corresponding  constants  are  identical.  The  assumption  that  »■(•)  and 
\>(.\  are  constant  functions  with  values  p and  x-  yield  the  representation 

C(nQ,n)  n( l -e  ) v n t N 

let  xp  be  any  process,  lhen  the  restriction  z to  (d  ,d  ♦$)  of 
the  frozen  bonding  continuation  of  y up  to  d *6  (if  it  exists) 
must  be  the  unixiue  function  satisfying 
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i = C(n(d<p).  z-n(d^));  z^)  = 

on  [d  ,d  +6)  or,  equivalently,  be  such  that  r,  = z - n(d  ) is  the 

iP  ifi  \p' 

unique  function  satisfying 

C = u(l-euC);  C(dv)  = Vldj-nidj  (4.2) 

on  [d^,d^+6).  The  smoothness  of  f.  - u(1-eu^)  suffices  to  ensure 
that  any  solution  to  (4,2)  on  [d^,d^+6)  is  necessarily  unique. 

Let  t _>  d^,  then  some  entirely  elementary  calculations  show  that 
the  function  e with  definition:  e(t)  = 0 if  C(d^)  = 0 and 

.“flO  . (l-e-^fV))"''  <"-3) 

if  C(d^)  / 0,  is  a solution  to  (4,2)  on  [d^,»).  Since  the  first  two 
terms  in  the  denominator  are  non  negative  and  the  remaining  term  is 
strictly  positive,  £(t)  is  well  defined.  Furthermore  when 
t(d^)  > 0 (resp.  t(d^)  < 0),the  replacement  of  by  unity 

shows  that  the  denominator  is  bounded  above  (resp.  below)  by  unity 
in  which  case  e(t)  +o(resp.  t 0)  as  t ®.  The  satisfaction  of 
AM^  is  immediate. 
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3.5  The  Evolution  of  Bonding  States 


In  order  to  characterize  the  bonding  (3,1)  corresponding  to 
each  process  ip  the  unique  function  satisfying  (3,2)  it  is  necessary 
for 

n = a{*  - M(l-eo(vi>-n  (5,1) 

subject  to  n(0)  = 1 , to  have  a unique  solution  n on  T0,d^)  or 
equivalently,  upon  setting  c * «p  - n t for 

L = <i>(l-a)  - ap(euC-l)  (5,2) 

subject  to  c(0)=0,  to  have  a unique  solution  t on  [0,d^). 

Let  <p  be  pieceuiee  linear  and  I = [i,f]  be  any  subinterval  of 
[0,d^)  on  which  <p  is  linear.  Then  (5,2)  has  the  solution 
L(t)  = L(i)  on  I whenever  £(i)  = 0.  Otherwise  (5,2)  integrates  to 
yield 

euC(t)  = (auu(t-i)  + e"°c(l)}  1 > 0;  tel  (5,3) 

whenever  c =(<J>(l-a)  + ap)/ap 

euC(t)  = c/{l-e‘Y(t'l)(l-ce'uf:(l))>;  t€l  (5,4) 

whenever  c / 0,  where  y - auuc.  Clearly  (5,4)  is  equal  to  (5,3)  in 
the  limit  as  c -*  0.  Since  the  first  two  terms  in  the  denominator  in 
(5,4)  are  non-negative  (resp  non-positive)  and  the  remaining  term  is 
strictly  positive  (resp  negative)  when  c > 0 (resp  c < 0),(5,4)  is 
well  behaved.  Also,  when  c > 0,  the  bonding  ultimately  lags  behind 
the  corresponding  stretch  by  !he  amount, 

C * Umt^uf.(t)  = (Ln  c)/o  (5,5) 

Since  the  right  hand  side 
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(t.O  e I *R  - ii>(t)  ( 1 -a ) + au(eu*-l) 

is  continuous  ami  locally  Lipschitz  in  L,  we  conclude  that  the 
solution  c corresponding  to  the  initial  data  L(i)  is,  in  fact, 
unique.  Iterating  this  result  we  have  the  desired  existence  and 
uniqueness  of  solution  to  (5,2)  for  all  piecewise  linear  processes. 

Let  \fi  be  the  continuously  differentiable  on  I = [ i , f ] . 

Clearly,  then,  we  can  construct  a sequence  of  piecewise  linear 

processes  ij>  such  that  ij>  -*  uniformly  on  I.  It  follows  that 
n n 

$ (t)  ( 1 -o)  + aVl(eot~l)  - i*>(t)(l-«)+«u(euC-l) 

n 

uniformly  on  I x R.  Since  for  each  n,  the  solution  r,  to 
in  - *n(t)(l-a)  + ap(euC-l);  tn(i)  * t(i) 

exists  and  is  unique  on  I,  it  follows  [2;  theorem  2.4]  that  there 

exists  a subsequence  r,  C such  that 

n ( 1 ) n ( 2 ) 

t(t)  = Lim  C ,Lv(t)  V t € 1 

k-*»  ,n*' 

is  a solution  to  (5,2)  on  I corresponding  to  the  initial  data  t(i). 
In  view  of  the  smoothness  of  the  right  hand  side,  we  conclude  that 
t is  unique.  Since  processes  are  continuous  and  piecewise  con- 
tinuously differentiable,  an  iteration  of  the  proceeding  result 
suffices  to  establish  that  (5,2)  has  a unique  solution  corresponding 
to  cacti  process.  The  characterization  of  the  evolution  of  bonding 
states  is,  then,  complete. 

Let  be  any  process  in  which  case  qi^  is  constant  on  [d  ,t] 

for  all  d v t < d «o.  Let  i.  be  the  corresponding  solution  to 
q>  q> 

ib, 2)  on  [d  , t)  . Since  the  value  of  c corresponding  to  <t>  = 0 is 

1 vp 

* 

unity,  we  have  that  f.  is  zero  and  hence  the  implication  : the 

bonding  corresponding  to  ultimately  tends  to  q>(d^)  from  above 

when  f. ( i ) < 0 and  helow  when  t(i)  > 0.  The  satisfaction  of  AM^ 

part  (a)  is,  then,  immediate. 
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3.6  The  Basic  Equations:  Reproduction  of  Experimental  Data 

We  assume  that  propellants  are  fading  memory  materials  in  tne 
sense  of  [ 1 ; sec  8].  Following  [3;sec  4]  and  fl;sec  l,l  together  with 
some  obvious  linearizations,  we  have 

ft 

S(*t)  - Ejn(g>  )-U  + E Mt)-n(ipt)]  + E G(t-i)$(i)di  + An(<p  ). 

o 

The  rapidity  with  which  a zero  stress  state  is  re-established  upon 
unloading  requires  that  =*  E$.  While  the  process  to  be  considered 
is  extremely  general  it  avoids  introducing  excessive  stress  relaxations 
in  which  case  the  stress  response  ought  not  to  be  critically  affected  by 
the  particular  choice  of  G(-)*  Eor  convenience  we  choose  G(t)  = tn 
where  -1  < n < 0.  The  explicit  governing  equations  are,  then, 

ft 

s(iPt)  3 E0  I ^(O  d'  + B0  "(lPt)  (M) 

and 

- «f*2>(t)-n(l-eu(<p(t)“n('pt))  I*.  n(4*0)  = 1 (6.2) 

Using  the  testing  approach  discussed  in  Section  2.0,  a propellant 
specimen  of  the  configuration  shown  in  Figure  1(a)  was  subjected  to 
the  characterization  test  described  in  Section  2.0.  The  series  of 
increasing  amplitude  sawtooth  events  (see  Figure  1(b)  at  a strain 
rate  of  0.0058  inches  per  second  was  followed  by  a sawtooth  event 
terminating  at  a strain  below  the  previous  maximum,  then  bv  sawtooth 
events  at  both  higher  and  lower  rates  and  to  various  amplitudes.  The 
actual  piecewise  linear  process  is  represented  in  Table  1.  Oata 
points  were  taken  from  the  tensile  tester's  pen  recorder  ("strip 
chart")  plots  of  load  cell  output  versus  time.  The  stretch  values 
(g>  (t))  were  calculated  from  the  crosshead  rate  and  time  (chart 
distance)  values  using  the  known  sample  length. 
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TARLt  1 


TEST  DATA  AND  CALCULATION  RESULTS  FOR  TP-H1011  PROPELLANT 


* 

t(secs) 

4>{  t) 

n(t) 

4>(tl(  ins  sec) 

Computed 

s ('P)  ( P . S . I , 

Measured 
.)  s(ip)(P.S 

0.0 

1 .000 

1 . 0000 

0.0 

0.0 

0 . 0 

0.70 

1 . 004 

1. 00  02 

0 . 00 1 8 

4.47 

5.20 

1. 70 

1.010 

1 . 0005 

0 . CO 58 

7. 4 7 

10.5  0 

2.70 

1.01b 

1. 00  08 

0 . 0058 

1 0 . t)  4 

16.25 

1.  70 

1.021 

1 . 00  1 1 

0.4058 

14.64 

2 1.65 

4.40  ^ 
4.76(1) 
5.20^' 

1 .025 

1.00  18 

0.0058 

18.64 

26.00 

1 . 027 

1 . 0021 

0 .0058 

21.32 

28 . 1 S 

1 . 025 

1. 00  23 

-0 . 00  6 8 

12.5  1 

14.30 

S.  70 

1 . 0 22 

1 . 00  2 1 

- 0 . o 0 5 8 

4.  1 0 

4.40 

t>.  70 

1.016 

1.00  22 

-0.0058 

4.5 

4 .75 

7.  70 

1.010 

1 .0020 

-0 . 0058 

1 . 1 • 

1 .60 

0.70 

1 . 005 

1 . 00  17 

-0 .0058 

-1.71 

-0.80 

4.  52 

1 . 000 

1.00  15 

-0 .0058 

-3.41 

-3.8  0 

11.22 

1 . 000 

1. 00  15 

0 . 0 

-1.00 

-1.10 

16.  22 

1 . 000 

1.0015 

0.0 

-0.  4 > 

-0.5  0 

21.22 

1 . 000 

1. 00  15 

0.0 

-0.24 

-0.30 

80  3.  22 

1.  0 00 

1 . 0006 

0 . i 

-0.0  1 

0.0 

004 . 12 

1.005 

1. 00  04 

0 .005 8 

5.06 

5.0  0 

80S.  12 

1.011 

1.00  12 

0.0058 

3.01 

'•.50 

806. 12 

1.017 

1. 00  15 

0.0058 

11.22 

12.50 

807.  12 

1.022 

1.0020 

0.C058 

15.34 

1 - . ’ 5 

808.  12 

1 .028 

1. 0028 

0 . 00 5 8 

21.80 

2 8 . 0 0 

804.  12 

804. 7S(  2 ) 

8 10.  12V"' 

1.0  34 

1.0044 

0.0058 

12.2  3 

34 .50 

1. 0 38 

1 . OOoJ 

0 . 0 0 5 8 

4 0.60 

33.  15 

1 . 0 16 

I.OOb-’ 

“0 .0053 

25.6  1 

21.  7 5 

811.12 

1.030 

1.00  74 

-0 . 0058 

1 1 . 2 7 

11.00 

3 12.  12 

1 .024 

1 . 00  74 

- 0 * 0 05 8 

7.77 

6.65 

813.12 

1.013 

1. 0072 

-0.0058 

4.1. 

1.3  5 

8 14.12 

1.012 

1 . OC  64 

-C . 0053 

1.14 

1.50 

8 1 5 . 1 2 

1 . 007 

1 . 00b7 

-0.0058 

-1.5- 

-0.65 

81o. 28 

1 . 000 

1 .006  1 

- 0 .0058 

- 4 .55 

- 1 . 1 •< 

318.25 

1.000 

1. 00b  1 

0.0 

-1.44 

-1.30 

8 15.  25 

1.000 

1.00  6 3 

0.0 

-0.40 

-0.2  5 

1375.25 

1.000 

1.0041 

0.0 

-0.04 

0.0 

1376.  15 

1.005 

1.00  44 

0 .0058 

4.44 

4 . -5 

1377.  15 

1.011 

1. 0047 

0.0058 

7.  8 7 

3.10 

1 178.  15 

1.017 

1 .00  50 

0.0058 

10.77 

11.50 

1 174.  15 

1.022 

1.0053 

0 . 0058 

14.13 

15.0  0 

1 380.  15 

1. 023 

1 . C054 

0 . 0058 

18.85 

14.75 

1 181 . 15 

1.034 

1. 00  b» 

0 . 0058 

26.6  0 

2'  .50 

1382.  15  _ 

1 383. 00 ( 3 ) 

1. 040 

1.0084 

0.0053 

18.4  1 

34. 2 > 

1.045 

1. 01  17 

0.0068 

4 4.?; 

45.00 

* Symbols  in  circles  identify  points  of  peak  strain 
(See  Figure  2). 
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TABLE  1 (CONT) 

TEST  DATA  AND  CALCULATION  RESULTS  FOR  TP-H1011  PROPELLANT 


1 303.  b5 

1.041 

1.01 30 

-0.0058 

25.41 

20.5 1 ' 

1384.15 

1.038 

1.0134 

-0.0058 

18.76 

15.50 

1305.  15 

1. 032 

1.01 35 

-0 .6058 

11.7t, 

1.65 

1 306. 15 

1.027 

1 . 0 1 34 

-0.0058 

7.66 

6.4  0 

1 383.  15 

1.015 

1.01 29 

-0 .1058 

1.68 

1.85 

1 30  9.  15 

1.009 

1.01  26 

-0.0058 

- 0 . 9 r- 

r .0 

1390. 75 

1 . 000 

1.01 22 

-0.0058 

-4.97 

-1.20 

1392.20 

1.000 

1.0122 

0.0 

-1.94 

-1.46 

1426. 20 

1 . 000 

1.01 20 

0.0 

-0.  3 3 

-0.2  0 

2026.20 

1.000 

1.00  88 

0.0 

-0.06 

0.0 

2028.00 

1.010 

1 .00  92 

0.0058 

7.4  9 

7.30 

2030.00 

1.022 

1. 0098 

0.0058 

11.06 

11.40 

20  32.00 

1.033 

1.01 08 

0.G058 

21.2  1 

2.7.5  0 

2033.50 

1.042 

1.0128 

0.0058 

34.74 

34.10 

2034.  50 
2035.00 ( 4 ) 
20  36.  50  W 

1 . 048 

1.0155 

0.0058 

47.9  1 

4 7.50 

1.055 

1.02  11 

0.0068 

61.99 

55.6  0 

1.051 

1.0228 

-0.0058 

30.7  0 

25.70 

2037.50 

1.046 

1. 0234 

-0.0058 

18.94 

15.50 

20  39.  50 

1 . 034 

1.02  32 

-0.C058 

9.94 

7 , / j 

2041. 50 

1.022 

1.0227 

-0.C058 

4.20 

3.4  0 

2043.50 

1.011 

1.0221 

-0.0058 

-0.9  2 

-0.10 

2045. 39 

1.000 

1. 02  16 

-0.0058 

-5.58 

-3.60 

20t>7 . 07 

1.000 

1.0215 

0.0 

-0.64 

0.30 

2697.07 

1.000 

1.0176 

0.0 

-0. 00 

0.0 

2698.57 

1.009 

1.0180 

0.0058 

6.61 

6.20 

2700. 57 

1.020 

1.0185 

0.C068 

11.89 

11.70 

2702.57 

1.032 

1.0191 

0.0058 

17.  2 3 

17.80 

2704.  57 

1 . 04  3 

1.0201 

0.0058 

25.53 

27.00 

2706.07 

1.052 

1.0222 

0.0058 

39.  54 

41.60 

2707. 07 
2708.47 

1.058 

1.0251 

0.0058 

52.7  2 

56.50 

1.066 

1.0313 

0.0058 

67.0  1 

65.0  0 

5709. "07 

1.062 

1.0329 

-0.0068 

37.0b 

3 3.20 

2710. 07 

1.057 

1.03  37 

-0.0058 

2 3 . 7 C 

19.8  0 

2712. 07 

1.045 

1.03  35 

-0 . 0G58 

14.14 

10.60 

2715.07 

1.028 

1.0327 

-0.0058 

5.  b4 

4.20 

2717.07 

1.016 

1.0322 

-0.0058 

0.60 

1.00 

2719.87 

1.000 

1.03  14 

-0.0058 

-6.  15 

-3.60 

2745. 22 

1.000 

1.0312 

0.0 

-0.8  1 

-0. 10 

3405.22 

1.000 

1.0271 

0.0 

-0.08 

0.0 

3407. 22 

1.012 

1.0276 

0.0058 

7.96 

7 .bO 

3409.22 

1.023 

1.0282 

0.0068 

13.06 

12.80 

3412.22 

1.040 

1.02  90 

0.0058 

20.  6 3 

21.80 

3414.22 

1.052 

1.0299 

0.0058 

28.  1C 

n .50 

3416. 22 

1.063 

1.0328 

0 . 0058 

46.86 

54.70 

3417.72  ^ 
3418.  38  (V) 

1.072 

1.0383 

0.0058 

o 5 . 4 3 

70.0  0 

1.076 

1.04  14 

0.C058 

71.0  1 

73.70 
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TABL  L 1 (CONT) 

TEST  DATA  ANl)  CALCULATION  RESULTS  FOR  TP-H1011  PROPELLANT 


34  14.22 

1.071 

1. 04  34 

-o.oosa 

lb.  1 7 

13.80 

3420. 22 

1 . ObS 

1 .04  38 

-0  . C0‘>8 

25.7  4 

21.70 

3422.22 

1 .054 

1 . 04  3b 

-0  . OO'-S 

17.04 

12.70 

3425. 22 

1 . 0 3b 

1.0428 

-0 .0053 

8.64 

b . 0 0 

1420.22 

1.014 

1. 04  20 

-0 . 0058 

1.19 

1.30 

)4  1 1 . 54 

1.  000 

1. 04  10 

-0.0058 

- b . 7 0 

-3.60 

14  SO . 2 4 

1.000 

1. 040  » 

0.0 

-1.00 

-0.20 

1448.  24 

1 . 000 

1.03  75 

o.c 

-0.  1 1 

0.0 

1444, 34 

1.006 

1.03  78 

0 . 0058 

5.4  7 

4.30 

400  1 . 34 

i.oio 

1.0381 

0.0048 

10.77 

10.00 

4 0 V 4 . 34 

1.0  IS 

1.0341 

0 .00  58 

18.0  >■ 

r . 5 o 

40  0«>.  14 

1.047 

1.0347 

0 .C058 

2 2.4  0 

2 3.3  0 

4000 . 34 

1 . 050 

1 . 04  04 

0 .0048 

28.  7 1 

H.00 

40b  4 . 14 
4010.24 (7) 
40 10.04 

1 . 0b4 

1 .04  10 

0.0053 

3 3.55 

1 >.0b 

1 . Ob4 

1.04  14 

0 .0058 

40.  1 7 

4 1.50 

1 . ObS 

1.0422 

-0.0058 

28.2  2 

'9.30 

4011. 04 

1 . 0S4 

1.04  2 3 

-0 . 0048 

22. 0 m 

14.00 

40  1 1. 04 

1 . 040 

1 . 04  18 

-0.0048 

14.8b 

11.00 

4 0 1 b . 8 4 

1.031 

1. 04  10 

- 0 . 0 0 4 El 

t)  , 0 i 

4.80 

4014.04 

1.013 

1 . 04C2 

- 0 . C 0 4 8 

-0.84 

0.2  0 

4022. 14 

1 . 000 

1. 0 3 45 

-0 .0048 

-b.  3 t 

- 3 . b 0 

4040.  14 

1.000 

1.0344 

0.0 

-0.‘K, 

-0 . 3 0 

4b48  . 14 

1. 000 

1.0  3 5b 

0.0 

-0.  1 1 

0.0 

4b  >0 . 14 

1. 00b 

1.0354 

0.0024 

3.4  0 
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4 b S 2 . 1 4 

1.012 
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1. 02  I 

1 .0367 

0.0024 

11.14 
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1.035 
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0.0024 
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1 . 04b 
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4b  b 3.  14 
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0 . 00 24 

2b.  0 7 
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4 b 7 0 . 14 

4H  7 2.04(8^ 
4b 7 3.  10  v-' 

1 . Ob  3 

1.03  48 

0.0024 

10. 8 0 

34.14 

1.0b4 

1.04  17 

0 .0024 

17.4  0 

4 4.44 

1 . Obb 

1.0424 

-0 .0024 

29.04 

2 7.74 

467b.  10 
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1 . 04  28 

-0.0024 
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18.74 
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1 . 048 
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4 b 0 1 . 10 
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-0 . 0024 
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1.0414 
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1.014 
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1.0403 

-0.0024 
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1. 0401 
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0.0 
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0.0 
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1. 047 
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3 3.28 
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4 0 . 5 1 
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TABLE 

1 (CONT) 

TEST 

DATA  AND 

CALCULATION 

RESULTS  FOR 

TP-H1011  PROPELLANT 

5324.14 

1.063 

1.0397 

-0.0115 

27.  3 C 

2 >.b5 

5 324.59 

1.057 

1.03  96 

-0.0115 

21  .3.’ 

19.50 

5325.59 

1.046 

1.0  3 92 

-0.0115 

13.42 

13.50 

532b. 59 

1.034 

1.0386 

-0.01 15 

7.34 

5.70 

5327.59 

1.023 
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1.0322 

0.0231 

14.44 
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1.0  22 

1.0327 

0.0231 

20.  38 

16.00 

6572. 22 

1 . 045 

1.03  38 

0 .0231 

31.82 

29.5  0 

6572.72 

1.057 

1.0344 

0.0231 

40.  0 7 

42.00 

657  3.27  (111 

1.069 

1.0359 

0.0231 

71.85 

o 3 . 2 0 

6573.71 

1.059 

1.0361 

-0.0231 
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1.0324 
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0.0006 
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17.30 
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1.0403 
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23. 9C 

25.50 
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1.051 

1.0421 
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15.40 

12.90 
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1.04  15 

-0.0006 

8.  32 
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1.000 

1.0395 

-0 . 0006 

-3.  1 5 

-1.50 

THE  MEAN  SQUARE  ERHOR  FOR  179  POINTS  IS  2.700 


The  force- responses 


tt  [o,  )Ks(  ) 

in  the  sixth  column  (Measured  s(g>)  ) were  calculated  from  the  chart 
readings  using  the  appropriate  factors  for  calibration  and  specimen 
cross-sectional  area. 

Each  trial  set  of  values^  0,  u > 0.  o 0,  and  n-  0 determines 
corresponding  values  for  e and  R0  from  the  requirement  that  the  root 
mean  square  error  of  the  force  responses  as  computed  from  (6,1)  and 
(6,2)  minus  the  corresponding  experimentally  measured  values  be  a 
minimum.  The  Bothm-Powel 1 functional  optimization  technique  [4] 
established  that  the  values 

-3.05 

a e 

n = 0. 132  x 10~3 

o=  )Q5 
n -0.103 

and  the  corresponding  values 

E * 495 

o 

Rp  7679 

rendered  this  minimum  equal  to  2.7  pounds  per  square  inch.  Clearly, 
the  tabulated  values  of  ,n(ip^)  in  Table  1 are  entirely  compatible  with 
our  intuitive  expectations  (also  see  the  "hysteresis"  plots  of 
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g'(t)  in  figure  i’).  In  coinin',  urn  >1  must  In  ompha.i.od 


(hat  this  study  is  one  of  an  interim  natunv  I ho  ina<  i in  ,n  ies  mi 

the  force- response  reproduction  appear  traceable  to  (he  usage  o( 
the  power  law  . This  usage  was  justly  critui/ed  Pv  uirtin  in  l *■  1 
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NOTE:  Symbols  in  circles 
denote  the  relative  time 
corresponding  to  peak  strain 
See  Table  1 for  actual 
elapsed  time. 


— I 1-  i r T 1 T 1 I T 1 I 

00  0 TO  0.4(1  0 SO  0 Bo  0.70  0.00 

STRAIN  ( S f R L T C H - I)  * 1 0 ' 


FIGURE  2.  STRAIN  AND  GONOING  HISTORY 


4.0  CONCLUSIONS 

The  work  represented  in  this  report  is  considered  to  he  a start 
ing  point,  with  a very  firm  foot i mi  in  continuum  mochaniis,  lor  a 
usable  constitutive  theory  for  solid  propellants.  A'  indicated  in 
Section  3.0,  the  deviation  of  calculated  force-response  (or  stress) 
from  the  measured  values  is  sufficient  to  be  less  than  totally 
satisfactory.  Dr.  Quinlan  suspects  the  power  law  represent  at  ion  used 
for  the  linearly  viscoelastic  component  of  the  force  response  as 
being  the  cause  of  this  deviation.  However,  a factor  which  cannot 
be  ignored  is  the  non-continuous  nature  of  solid  propellants,  whiih 
may  require  introduction  of  non-continuum  considerations  before  a 
total ly  workable  propellant  constitutive  law  can  he  produced 

The  degree  to  which  the  results  herein  reproduce  damaged  propel 
l ant  response  is  extremely  encouraging,  and  further  work  in  tins 
area  will  undoubtedly  resolve  the  remaining  problems. 
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APPENDIX  - MATHEMATICAL  NOTATION 


Section  3 of  this  report  uses  certain  mathematical  notation  ami 
terminology  which  may  be  unfamiliar  to  many  reader',  of  this  report. 

This  appendix  is  intended  to  aid  these  readers  in  following  the 
theoretical  development.  With  the  explanations  given  in  this  appendix, 
the  remaining  development  of  terminology  and  use  of  symbols  in 
Section  3 should  be  self-explanatory. 

The  notation  used  is  explained  below  using  a tabular  format,  in 
no  particular  order.  Notation  which  should  be  clear  to  most  engineers 
is  (hopefully)  not  included. 

NOTATION  MEANING 

: "Such  that";  e.g.,  t:(t-l)=0  means 

"t  such  that  (t-1)  = 0". 

e "Is  an  element  of  (a  set)";  e.g., 

a e (A)  means  "a  is  an  element  of  the 
set  (A)".  Also  can  be  read  "belongs  to". 

V "For  all". 


e "Is  a subset  of";  e.g.  C c p means 

" C is  a subset  of  P" 


NOT A I ION 

4 or  t 

-* 

or 


La,b] 

La.b) 

(a,b] 

(a.b) 

(•) 

(A.B.-) 


MEANING 

"Approaches  from  above"  or  "approaches  from 
below";  e.g.  , t 4 0 means  the  same  as 
"t  - 0;  t > 0". 

"Maps  into";  e.g.,  f : Uf  -►  Rf  means 
"f  such  that  Of  maps  into  Rf"  where  Of  is  the 
domain  of  the  function  f and  Rf  is  the  range 
of  values  assumed  by  the  function.  Also, 
x m f(x)  can  be  read,  "Assign  to 

x (in  some  obvious  domain)  a function  f(x)"; 

2 

e.g.,  x x is  more  fully  written  as 

( x ) 2 : (-“■',  co  ) _ ] 

A closed  interval , which  includes  a and  b ; 
e.g.,  x : [a,b]  means  the  same  as  x:  a < x > b. 

An  interval  which  includes  a but  excludes  b. 

An  interval  which  includes  b but  excludes  a. 

An  open  interval,  which  excludes  a and  b. 

Means  "A  variable  which  is  obvious  from  the 
context". 

Means  the  variable  represented  by  (•)  is  allowed 
to  vary  while  A and  15  are  held  constant. 
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NOTATION 


MEANING 


! 


resp. 


" Respectively" . 


sup  Supremum,  or  least  upper  bound. 

moduluo  "Except  for". 

Lipschitz  A function  which  is  "Lipschitz"  has  a continuous 

derivative.  "Locally  Lipschitz"  is  a more 
stringent  restriction  than  "continuous"  and 
a less  stringent  restriction  than  "continuously 
differentiable". 


AMp  AM£,  NRF  Reference  to  an  axiom  (underscored  in  the 

statement  of  the  axiom)  or  an  assumption 
(underscored  in  the  statement  of  the 
assumpti on) . 


I x R (As  used  in  the  equation  following 

Equation  (5,  5)  , i.e.,  (t,  I )£  I x R ) 

an  interval  of  t (defined  previously  as  a 
subinterval  of  [0,d^]  ) and  a corresponding 
interval  of  5 consisting  of  all  real  numbers. 
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